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Successful realization of spin-orbit coupling in atomic gases by the NIST scheme opens the prospect of 
studying the effects of spin-orbit coupling on many-body physics in an unprecedentedly controllable way. Here 
we derive the linearized hydrodynamic equations for the normal atomic gases of the spin-orbit coupling by the 
NIST scheme with zero detuning. We show that the hydrodynamics of the system crucially depends on the 
momentum susceptibilities which can be modified by the spin-orbit coupling. We reveal the effects of the spin- 
orbit coupling on the sound velocities and the dipole mode frequency of the gases by applying our formalism to 
the ideal Fermi gas. We also discuss the generalization of our results to other situations. 


The persisting quest to simulate charged particles in solid state systems by neutral atoms'^^ stimulated the pioneering ex¬ 
perimental achievement of realizing “spin-orbit” coupling in atomic gases through the Raman process by the NIST group"^. In 
the NIST scheme, a small magnetic field in the z direction was used to open the degeneracy of atomic hyperfine spins, and two 
Raman lasers aligning in the x direction shed on a gas of Bose atoms gave rise to a coupling bilinear in momentum and (pseudo-) 
spin of atoms in one direction. After a unitary transformation^, the resulting single atom Hamiltonian has the form (we take 
h = 1 throughout) 
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where m is the atomic mass, cti are the atomic (pseudo-) spin operators, is the Raman laser wave vector, S is the detuning 
tunable by changing the frequency difference between the two Raman lasers, and Cl is the Rabi frequency for the Raman process. 
Note that the spin-obit coupling in Eq. (1) is different from the Rashba and the Dresselhaus forms; when the lasers are turned 
off, i.e., H = 0, the coupling kxUz in the kinetic energy can be eliminated by a gauge transformation^. Only nonzero Cl gives 
rise to nontrivial spin-orbit coupling. Experimentally k^ is fixed by the wavevector of the Raman lasers and Cl can be tuned by 
the laser intensities. Direct diagonalization of Hq yields the single atom dispersions 



with + (—) standing for the upper (lower) branch. The same scheme was later applied to Eermi gases®’^. 

The ground state of a Bose gas subject to the spin-orbit coupling realized by the NIST scheme can be either the stripe, the 
magnetic, or the non-magnetic states depending on the spin-orbit couplings and interatomic interactions'*’^’^. Recently the finite 
temperature phase diagram of the Bose gas was determined experimentally for zero detuning b = 0®. A subsequent perturbative 
calculation reproduced the correct trend of the thermal effects on the phase boundary between the stripe and the magnetic 
phases 

The effects of the spin-orbit coupling of the NIST scheme on the dynamics of atomic gases were first experimentally investi¬ 
gated through the collective dipole oscillation of a Bose-Einstein condensate in a harmonic trap* * . In the absence of the spin-orbit 
coupling, the atomic gas has the Galilean invariance, which guarantees the dipole oscillation frequency Wd equal to the harmonic 
trapping frequency The dispersions (2) given rise to by the NIST scheme apparently break down the Galilean invariance, 
which indicates that ujd can be different from uiq. Nevertheless, since during the small dipole oscillation the Bose-Einstein con¬ 
densate accesses mostly the lower branch states around a minimum of e_ k at momentum kg, it is sufficient to approximate the 
dispersion e_.k ~ e-.ko + (k —ko)^/2m*, whose form restores the Galilean invariance though the “effective mass” m* depends 
on the spin-orbit coupling. The experimental data for the dipole frequency of the Bose-Einstein condensate with the spin-orbit 
coupling can be mainly explained by the “effective mass” approximation Wd ujq a/ m/m*' ' . 

Up to now, the study of the collective modes and the hydrodynamics of the spin-obit coupled atomic gases mainly focuses on 
the Bose-Einstein condensates, in which case the existence of a single condensate wave function greatly simplifies the theoretical 
treatment*’* However, to describe the dynamics of Bose gases with a substantial normal fraction and Eermi gases in the 
presence of spin-orbit coupling requires a more general framework. In this work, we consider atomic gases in the hydrodynamic 
regime and derive the linearized hydrodynamic equations for the normal atomic gases with the NIST spin-orbit coupling for zero 
detuning (5 = 0. We show that in the absence of the Galilean invariance, the hydrodynamics of the system crucially depends on 
the momentum susceptibilities, which can be modified by the spin-orbit coupling. We apply our general formalism to the ideal 
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Fermi gases and reveal the effects of the spin-orbit coupling on the sound velocities and the dipole oscillation frequency of the 
gases. 

Results 

Linearized Hydrodynamic Equations. The normal atomic gases with the NIST spin-orbit coupling for zero detuning b = 0 
conserve the number of atoms N, the energy E, and the (pseudo-) momentum K [cf. Eq. (1)]. The three conservation laws in 
their differential forms are 




^ + v.n = o, 


(3) 
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Here e, g and n are the local densities of energy, momentum and number of atoms respectively, and IT and j" are the energy 
current, momentum current tensor and number current respectively. We assume that the inter-particle collision is so frequent 
that the system is in the hydrodynamic regime. The density of any physical quantity O can be calculated by the local grand 
canonical ensemble with the distribution exp[—/3(iF — piV — K • v)], where /3, p and v are the inverse of temperature T (we take 
fee = 1 throughout), chemical potential and velocity fields which are functions of position and time. The total Hamiltonian H 
includes Hq and interatomic interactions. Note that due to the spin-orbit coupling, the spin degrees of freedom are not conserved. 
Spin dynamics in the presence of spin-orbit coupling has been studied in the context of both electron gases and atomic Fermi 
gases 

We focus on the dissipationless limit. The condition that the entropy change of the total system is zero determines the 
constitutive relations ^^42 


r = «v (6) 

Hab = P^ab “b '^a9b (7) 

= (e + P)v, (8) 

where P is the pressure and the subscript stands for the component index. If the atoms are subject to an additional external 
potential C/ext (r), the momentum conservation equation becomes 

^+^V„(P<5a6+t;a56)+^Vf,C/ext = 0. (9) 

a 

The absence of the Galilean invariance manifest through Eq. (2) is a key feature of the atomic gas with spin-orbit coupling 
engineered by the Raman processes driven by external lasers. For normal atomic gases with the Galilean invariance, the momen¬ 
tum density g is always equal to mnv for arbitrary velocity field v. However, this equality generally breaks down in the presence 
of spin-orbit coupling. The consequence of the lack of the Galilean invariance has been shown to affect the superfluidity and the 
bright solitons in Bose-Einstein condensates^^44 Nevertheless, if we are interested in small variations of /3, p and v away from 
the global equilibrium, we can keep to the first order of the variations and have 


5a(/3,/i,v) =XahU6, (10) 

e(/3,M,v) =i^Xat,UaUf, + e(/3,/r,v = 0), (11) 

a,& 

where the momentum susceptibility is 

Xab = 5ga/5vb\^=o- (12) 

For simplicity, we have assumed that g(/3,/ 2 ,v = 0) = 0, i.e., when the gas is not moving, its momentum is zero. This 
assumption is true for the NIST scheme with zero detuning. Generalization of our results to cases without this assumption is 
straightforward. Apparently Xab must be symmetric. We choose the coordinate system such that Xab = XaSab- Using Eqs. (10) 
and (11), we distill the linearized dissipationless hydrodynamic equations into 


d'^Sn 
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where 5n is the variation of number density n away from its global equilibrium value no and {dP/dn)s is taken at fixed entropy 
per particle. At T = 0, Eq. (13) further simplifies into 


d'^Sn 




Xa 



(14) 


The phenomenological hydrodynamic equations (13) and (14) shall be applicable to both normal Bose and Fermi gases with the 
spin-obit coupling. 

Momentum Susceptibility. The linearized hydrodynamic equations for atomic gases with spin-orbit coupling depend explicitly 
on the momentum susceptibility Xab which can be calculated by the formula^^ 


{l\Kb\f){f\Ka\l) 

V Ef-Ei 
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(15) 


where Ei and |i) are the eigenvalues and eigenstates of the many-body Hamiltonian H, pi is the distribution function, and V is 
the gas volume. 

To manifest the effects of the spin-orbit coupling on Xab, we assume that the inter-atomic interaction is weak enough that we 
can evaluate Xab using the ideal gas Hamiltonian, and find 


Xab — 





^b^af (^a.k) : 


(16) 


where /(e^.k) is the Fermi (Bose) distribution function for fermionic (bosonic) atoms, D is the dimension of the gas. The 
information of the spin-orbit coupling is encoded in the single atom dispersion k- For degenerate Fermi gases, Xab shall be 
dominated by the contribution from close to the Fermi surface. 

We calculate Xxx (Xx) by Eq. (16) for ID, 2D and 3D Fermi gases with spin-orbit coupling generated in the x direction by 
the NIST scheme with <5 = 0. The (quasi-) ID or 2D gases can be achieved when there is strong confinement in the y direction 
or in both the y and z directions. For the ID Fermi gas of chemical potential p at zero temperature, we find 
1) pj £r ^ 1 
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(18) 


Xxx _ ^ + 1 + \/(/^/£r 1)^ (ll/2er)^ 

mn~ 2[p/e, + (H/4e,-)2] ’ ^ ^ 

where the recoil energy is e, = k//2m. 

Figure (1) shows at zero temperature how Xxx for the ID Fermi gas changes with H for different atomic densities n. The 
plateaus appearing at small H are peculiar to this ID case and can be understood in the following way. When <5 = 0 and H = 0, 
the two energy bands given by Eq. (2) touch at energy e, when kx = 0. The Fermi surface crosses the lower band at four points 
for n/k, < 2 /tt{k, 0.64), and crosses both bands each at two points for n/k, > 2/7r. When is increased from zero, an energy 
gap ~ H opens at kx = 0. One can show analytically that for small densities n before the point of the lower band at = 0 
becomes lower than the Fermi surface, or for large densities n before the upper band bottom at = 0 becomes higher than the 
Fermi surface, Xxx/mn =1-1- {2k,/ Tm/^. The further away the density n is from 2/tx, the larger H the plateau persists to. Of 
course for our hydrodynamic approach to be applicable, H must be sufficiently larger than local equilibration rates. In the large 
H limit, the Fermi surface crosses the lower band at two points and Xxx/mn —> (1 — Ae,/Vl.)~^ [cf. Eq. (2)]. 

We also plot Xxx versus atomic density n for the ID Fermi gas with different O in Fig. (2). When the density is high, the Fermi 
surface lies at high momenta where the spin-orbit coupling has negligible effects, and Xxx approaches mn as expected. When 
Vl/e, = 2, has minima at two distinct kx- For Vl/e, = 4, has zero curvature at its single minimum. These features 
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of e_^fc^ give rise to the corresponding divergences of Xxx/mn in the low density limit shown in Fig. (2). When O/cr = 5, the 
ratio Xxx/'mn, though finite, is enhanced to be substantially larger than unity by the spin-orbit coupling. Similar behavior of 
Xxx would be found in the 2D and 3D Fermi gases as well. 

In Fig. (3) is shown the finite temperature behavior of Xxx for the Fermi gases with various Rabi frequency 17 in different 
dimensions. The spin-orbit coupling moves Xxxl'mn more away from unity at lower temperatures or in lower dimensions. When 
T is finite, the low density limit corresponds to the chemical potential fx —^ —oo. In this limit, to zero order, the distribution 
function / in Eq. (16) can be approximated by the Boltzmann distribution function; Xxxl'mn acquires the same value for fixed 
17 and T in different dimensions. 

Sound Velocities. In the case that there are no external potentials, i.e.. C/ext = 0, we can read off the sound velocities along the 
principal axes fromEq. (13) as 


2 ( 9P \ 


( 20 ) 


where s is the entropy per particle. Besides the momentum susceptibility Xa, the sound velocities also depend on the equation 
of state via the adiabatic compressibility Kg = {dn/dP)^ jn. Eor convenience, we recast the sound velocities into 


1 Cp 
XaK-T Cv ’ 


( 21 ) 


by using the thermodynamic identity Ks/kt = CvjCp, where kt = jr? is the isothermal compressibility, 

Cy = {dE/dT)y jy and Cp = {dE/dT)p = Cy + VTup i9P/dT)y are the specific heats. 

We calculate Cx by Eq. (21) combined with our previous results for Xxx for ID, 2D and 3D ideal Eermi gases with spin-orbit 
coupling generated in the x direction by the NIST scheme with (5 = 0. Eor convenience, the sound velocity is normalized by 
Co which is the corresponding sound velocity of the ideal Eermi gas without the spin-orbit coupling at the same density. Eigure 
(4) shows that at finite temperatures Cx for the ID case changes non-monotonically as the density of the Eermi gas varies. This 
non-monotonic behavior can be understood from the zero temperature limit. At zero temperature, according to the Gibbs-Duhem 
relation we have dP/dn = ndji/dn', from Eq. (20) the sound velocities depend on the density of states at the energy scale of 
the chemical potential ji. In ID, when the chemical potential approaches the minimum of the NIST dispersion e+,k. the density 
of states is divergent. This divergence suppresses the sound velocity to zero as shown in Eig. (4) for f7/er = 4, and gives rise to 
a discontinuity there. Therefore at finite temperatures, the thermal effects smear out the singular behavior of at T = 0 and 
result in the non-monotonic one shown in Eig. (4). When the dimension of the gas is raised up to 2D and 3D, the effects of the 
density of states on Cx persist and cause a finite jump in 2D shown in Eig. (5) and a cusp in 3D as in Eig. (6) at zero temperature. 
Thus the finite temperature behavior of Cx follows the general trend of its zero temperature one. 

Dipole Mode in Harmonic Traps. Collective modes of atomic gases confined in harmonic traps ?7ext = ma;Qr^/2 are important 
physical observables. To investigate how the spin-orbit coupling affects the dipole mode frequency Wd of the normal atomic 
gases with the NIST spin-orbit coupling, instead of solving Eq. (13), we adopt a variational formalism which is equivalent to the 
linearized hydrodynamic equation (13)^^“^^. We start with the action A = J dtL and the Lagrangian is 


L= J d^r ^ Xavl nUe^t + < 

K a 

ds 


dn ^ . 

^ + V-(nv) 


+r] 


dt 


-f V • (sv) 
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where s is the entropy density, and cj) and rj are the Lagrangian factors introduced to enforce the number conversation dn/dt + 
V • (nv) = 0 and the dissipationless condition ds/dt + 'S/ ■ (sv) = 

To estimate the frequency of the dipole mode oscillating in the principal axis x direction, we assume the ansatz n{r,t) = 
no(r — xo(t)x) and s(r, t) = so(r — xo{t)x) with uq and sq the functions at equilibrium. This ansatz is motivated by the fact 
that the dipole mode is mainly the “center of mass” motion of the gas. The number conservation gives v = dxo{t)x/dt, which 
also maintains the dissipationless condition. After substituting the above ansatz into the Lagrangian, we have 


L=Ln 



1 

2 



2 2 
mujQXo, 


(23) 


where Lq is the Lagrangian at equilibrium. Erom Eq. (23), we obtain 


d^xo m f d^rno 9 


( 24 ) 
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The dipole mode frequency is 


uja /rnf d^rno\^''^ 

ujo [ f d^rxx ) 


(25) 


The independence of Wd on the equation of state originates from the ansatz we use. It is worth mentioning that within the local 
density approximation, the result (25) agrees with the one given by the sum rule approach^^. Figure (7) plots the dipole mode 
frequency of ID, 2D and 3D normal Fermi gases of iVtot number of fermions with the spin-orbit coupling by the NIST scheme 
with (5 = 0. Due to the enhancement of Xx compared to mn shown above, the dipole mode frequency Wd/wo is generally 
suppressed below unity. 

Discussion 

The hydrodynamic equations (13) and (14) are valid for normal atomic gases with any spin orbit coupling under the condition 
that the momentum K = 0 if the velocity v = 0. Apparently even within the NIST scheme, if (5 7 ^ 0, this condition does not 
hold. Equation (10) shall include first order variation of p and /?, and Eq. (11) shall change correspondingly; the generalization 
to Eqs. (13) and (14) shall be straightforward. We have revealed the effects of spin-orbit coupling on the hydrodynamics of 
atomic gases by explicitly applying our general formalism to ideal Eermi gases. Interatomic interactions are expected to make 
quantitative changes to the results presented above. However, to take into account the interaction effects, one needs reliable 
determination of the susceptibilities Xab and the equation of state of the gases, which is beyond the scope of our work. The 
generalization of our hydrodynamic equations to the cases in which there is condensation requires correct treatment of the 
“superfluid” part. Since the form of the “superfluid” part depends on the exact structure of the order parameter, we leave the 
generalization to a future study. 

Our calculation manifests the effects of the spin-orbit coupling on physical observables such as sound velocities and dipole 
mode frequency. Previously the sound velocity of a unitary Eermi gas has been measured by creating a local density variation 
through a thin slice of green laser and monitoring the propagation of the density variation afterwards^®. Recent measurement of 
second sound velocity in the superfluid phase of unitary Eermi gases has also been achieved^®. The dipole oscillation of a Bose- 
Einstein condensate with spin-orbit coupling was excited by a sudden change of the Raman detuning and yielded clear violation 
of Kohn theorem". We expect that similar experimental techniques can be employed to confirm our theoretical predictions. 
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FIG. 1 : Momentum susceptibility of the ID ideal Fermi gas at zero temperature versus Rabi frequency fl. Dependence of the 
momentum susceptibility Xxx of the ID ideal Fermi gas on D for different gas densities n. 
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FIG. 2: Momentum susceptibility of th® ID ideal Fermi gas at zero temperature versus atomic density n. Dependence of the 
momentum susceptibility of th® ID ideal Fermi gas on the Rabi frequency fl and the gas density n. The grey dash line is x^^ = tnn for 
cases in which the spin-orbit coupling has negligible effects. 
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FIG. 3: Momentum susceptibility x^x of the ID, 2D and 3D ideal Fermi gas at finite temperature versus atomic density n. Dependence 
of Xxximn on the temperature T and the density n for the ID, 2D and 3D ideal Fermi gases. 
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FIG. 4: Sound velocity Cx of the ID ideal Fermi gas versus atomic density n. Dependence of sound velocity Cx of the ID ideal Fermi 
gas on the temperature T and the density n. The sound velocity is normalized by cq the corresponding sound velocity of ID ideal Fermi gas 
without the spin-orbit coupling at the same density. The green solid line is for T = 0, the blue short-dash line for T = O.Scr, the red dash line 
for T = l.Ofir, the purple short-dash-dot line for T = 2.0er all with D = 4er. The black solid line is for T = l.Oer with D = 5er. Inset is the 
density of states in unit of er at D = 4er. 
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FIG. 5: Sound velocity Cx of the 2D ideal Fermi gas versus atomic density n. Dependence of sound velocity Cx of the 2D ideal Fermi gas 
with D = 4er on the temperature T and the density n. Here co is the finite temperature sound velocity of the 2D ideal Fermi gas without the 
spin-orbit coupling. See the caption of Fig. 4 for the temperatures for each curves. Inset is the density of states in unit of mV 2 D at D = dcr- 
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FIG. 6: Sound velocity Cx of the 3D ideal Fermi gas versus atomic density n. Dependence of sound velocity Cx of the 3D ideal Fermi gas 
with D = 4er on the temperature T and the density n. Here co is the finite temperature sound velocity of the 3D ideal Fermi gas without the 
spin-orbit coupling. See the caption of Fig. 4 for the temperatures for each curves. Inset is the density of states in unit of mYsoki at D = de,. 
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FIG. 7: Dipole mode frequency uJd versus temperature T. Dependence of dipole mode frequency at finite temperature on the temperature 
and the Raman coupling for the ID, 2D and 3D systems. We have taken cuq = 27r x 164 Hz and Cr = Svr x 8340 Hz® and the total number of 
particles Nm = 10^, 10^, 10®, respectively. 










